Interaction of a discrete breather with a lattice junction 
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We study the scattering of a moving discrete breather (DB) on a junction in a Fermi-Pasta-Ulam (FPU) chain 
consisting of two segments with different masses of the particles. We consider four distinct cases: (i) a light- 
heavy (abrupt) junction in which the DB impinges on the junction from the segment with lighter mass, (ii) a 
heavy-light junction, (iii) an up mass-ramp in which the mass in the heavier segment increases continuously as 
one moves away from the junction point, and (iv) a down mass-ramp. Depending on the mass difference and 
DB characteristics (frequency and velocity), the DB can either reflect from, or transmit through, or get trapped 
at the junction or on the ramp. For the heavy-light junction, the DB can even split at the junction into a reflected 
and a transmitted DB. The latter is found to subsequently split into two or more DBs. For the down mass-ramp 
the DB gets accelerated in several stages, with accompanying radiation (phonons). These results are rationalized 
by calculating the Peierls-Nabarro barrier for the various cases. We also point out implications of our results in 
realistic situations such as electron-phonon coupled chains. 

PACS numbers: 63.20.Pw, 63.20.Ry, 87.10.+e, 66.90.+r 



I. INTRODUCTION 

Static discrete breathers (DB) are time-periodic, persistent, 
intrinsic localized exact modes in nonlinear lattices. Rigor- 
ous proofs of their existence have been obtained and system- 
atic studies of their properties were carried out using various 
(approximate) complementary approaches, see e.g. Ill for an 
overview. In contrast, as first noticed in numerical investi- 
gations and then justified theoretically, moving DBs exist as 
approximate solutions in nonlinear lattices, both Hamiltonian 
and non-Hamiltonian (with dissipation and periodic forcing). 
These solutions are known to be rather stable (i.e., long-lived) 
and have been an object of constant investigation during the 
last decade, see [|| for a non-exhaustive list. 

Different physical systems in which there are realizations 
of (moving) DBs include conjugated polymers ^[j}, charge- 
density-wave materials (e.g., metal-halogen electronic chains 
|^]), Josephson ladders [0], coupled electron-vibron lattice 
systems |Q| and spin chains m\. Sputtering on crystal surfaces 
and damage tracks in certain mica minerals have also been 
attributed to moving breathers Experimentally, breathers 
have been probed by ultrafast resonance Raman [^|] and in- 
elastic neutron scattering [ 1 1 among other techniques. 

In a recent series of papers [flip, the problem of the inter- 
action of a moving DB with an impurity was addressed in the 
case of a lattice with nonlinear on-site potential and harmonic 
first-neighbor coupling. As it was shown, this interaction can 
lead to reflection, transmission or trapping of the DB at the im- 
purity, depending on the initial velocity, amplitude and phase 
of the DB, as well as on the 'strength' and spatial extent of the 
impurity. 

Our objective here is to investigate the scattering of a DB at 
a junction in a (nonlinear) FPU chain consisting of two seg- 
ments that are "slightly different"-i.e., for instance, with dif- 
ferent interaction parameters or with different masses of the 
particles in the two segments. The reason for choosing the 
FPU chain is simple. It has historically provided a testbed 
for exploring novel nonlinear phenomena in discrete systems. 



In addition, it is one of the simplest nonlinear (polynomial) 
potentials amenable to some analytical calculations. 



Our preliminary numerical simulations indicate that these 
two types of problems are not qualitatively very different. 
Therefore, here we will concentrate exclusively on the second 
type of configuration, the one with slightly different masses 
on the two sides of the chain. A physical realization of this 
configuration could be in low-dimensional electronic materi- 
als with different electron-phonon coupling or two segments 
with different isotopes (e.g., carbon isotopes in conjugated 
polymers [[|^]] and platinum isotopes in metal-halogen chains 
[^]), Josephson- Junction arrays [g] with dissimilar interaction 
strengths, optical fibers with two different refractive indices 
Jl2|], etc. We note that the scattering of Toda solitons at a 
mass interface was studied previously [p"3|]. To the best of our 
knowledge, the scattering of a DB at such an interface has not 
yet been investigated. 



The paper is organized as follows. In Sec. II we present 
the details of the FPU model in a homogeneous chain, an es- 
timate of the Peierls-Nabarro barrier for moving DBs, and fi- 
nally, some details on the numerical initialization of a moving 
DB. Section III contains results for both the light-heavy and 
the heavy-light mass junction, where we elaborate on the re- 
flection and transmission (and eventually on the splitting) of 
the DB. Interaction of the DB with both the up mass-ramp and 
the down mass-ramp is discussed in Sec. IV, where we explore 
DB reflection (with eventual trapping) and acceleration (with 
eventual splitting). In Sec. V we summarize our main find- 
ings and enumerate some of the open questions. Details of the 
Peierls-Nabarro barrier calculation, using a new perturbative 
technique, for the various homogeneous and inhomogeneous 
cases are relegated to an Appendix. 
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II. THE MODEL 



A. The FPU model 



The FPU model represents a one-dimensional (ID) chain 
of particles with no on-site potential (i.e., an acoustic chain), 
with the Hamiltonian 
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where a and (3 denote, respectively, the strengths of the linear 
and nonlinear nearest-neighbor interactions; a represents the 
lattice constant (i.e., the equilibrium distance 
between neighboring sites), and m is the mass of the parti- 
cles. For simplicity, all these quantities (and those we will 
introduce later) are expressed in dimensionless units. The cor- 
responding equation of motion for a generic particle is: 

mx n = a (x n +i + x n -i - 2x n ) + ft [fa n - x n -\ - a) 3 
- (x n+ i -x n - a) 3 } . (2) 

In terms of the elongations u n = x n — x n q , it becomes: 

mu n = a (u n +l + U n -1 - 2u n ) + (3 [{u n - U n -i) 3 

- {Un+1 - U n ) 3 ] , (3) 

or, by introducing the relative elongations of neighboring sites 



(x n x n) fan— I x n-l) 



(4) 



m,T n 



a (r, 



n+l 



2t„ 



01 



'n+l 



2t 3 



(5) 



As it was shown (see, for example [|lj] for an overview, and 
references therein, and [|l4|]), the FPU lattice admits DB-like 
solutions (stationary, localized, time-periodic modes) with pe- 
riods Td b that are smaller than the minimum period of the 
phonon spectrum, i.e., 



Tdb < 7T y/m/a . 



(6) 



Also, as shown, for example in [|15| and Jlq], the most local- 
ized of these modes are an odd- type mode with an "approxi- 
mate" pattern of the amplitudes of the elongations u n of the 
form: A dd{---, 0, —1/2, 1, 1/2, 0, ...), and an even-type mode 
JT^ ] : A even (..., 0, — 1, 1, 0, ...). The amplitudes A are deter- 
mined by the interaction constants a and (3 and by DB's fre- 
quency ojdb = 2-k/Tdb- For given interaction constants, the 
A's decrease with increasing Tdb- On the contrary, keeping 
Tdb fixed and decreasing the interaction constants generally 
leads to an increase in the amplitudes A. By "approximate" 
above we mean that (as seen in flip]]) these patterns are exact 
only for a pure even-order anharmonic lattice in the limit of 
increasing order of anharmonicity. Nevertheless, only minor 



corrections are needed in order to make these patterns 'more 
precise' solutions of the FPU lattice, their symmetry being 
preserved. Mainly, these corrections refer to the fact that the 
DB can extend over more than three, and two sites, respec- 
tively, for odd and even modes. Although, in theory, a DB has 
an infinite extension [with an exponential decay of the ampli- 
tude of the relative elongation as one moves far away from 
the center (maximum amplitude sites) of the DB], in practice, 
however, one can restrict the analysis to five, and four sites, 
respectively, for the two types of modes of the DB mentioned 
above. 

To evaluate the relative elongations for the two configu- 
rations, and their corresponding energies, we introduced a 
simple perturbative technique that uses the ratio between the 
square of the maximum phonon frequency and the square of 
the DB frequency as the perturbation parameter: 
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combined with a rotating wave approximation, RWA (see, 
e.g., [[ij and references therein). The results of our calcula- 
tions, presented below and, in more detail, in the Appendix, 
can be compared with the numerical results of Green's func- 
tion method (that is also based on RWA). For example, for the 
even-symmetry mode, our calculations [up to 0(e 2 )] agree 
generally up to an error of no more than 4% with the results 
of [|l8]] obtained with Green's function method. The error in 
evaluating the configurations (as compared with the results 
of the exact numerical method of the analytical continuation 
from the anticontinuous limit [|l9|]) is essentially connected 
with the limitations of RWA, and therefore becomes progres- 
sively smaller for 'heavier' DBs, i.e., DBs that are progres- 
sively further away (in frequency) from the phonon band limit. 

The primary ingredients of the analytic method are the 
ansatz concerning the temporal evolution of particles' elon- 
gations: 



U n {t) = A^n COsfalDBt) ■ 



(8) 



(where A and are the amplitude and the shape func- 
tion, respectively), together with the RWA that entails ne- 
glecting higher-frequency harmonics [i.e., cos 3 (uiDBt) ~ 
(3/4) cos(woBi)]. Including these elements in Eq. (|j), one 
obtains an infinite set of nonlinear coupled equations for the 
shape function: 

mi^DB^n = «(2£„ - £ n+ i - £ n _i) 

+ - A PA 2 [{tn ~ in+lf + (& ~ £n-l) 3 ] . (9) 

Or, in terms of the relative elongations: 

r„(i) = A( n cos(w DB i) , ( n = £„ - f n _i , (10) 
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Next, we consider the following expansion of the shape func- 
tion in terms of the small parameter e, see Eq. (0): 



in = e 
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and then proceed through the usual steps of a perturbative cal- 
culation. For details of these calculations, refer to the Ap- 
pendix. 



1. The Peierls-Nabarro barrier for the homogeneous FPU chain 



As illustrated in Fig. 1(a) on an actual example, the DB 
translates from one lattice site to another by continuously de- 
forming its shape, alternately, between an odd-type of con- 
figuration and an even-type one. Therefore, in a discrete lat- 
tice there is an energy cost associated with moving a nonlin- 
ear localized mode by a lattice constant-this represents the 
so-called Peierls-Nabarro barrier (PNB), see p0[]. It can be 
estimated by calculating the energy difference between even- 
and odd-type configurations. The results presented in the Ap- 
pendix allow us to evaluate the PNB in an homogeneous chain 
(i.e., all particles with same a, j3, and m): 
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where the superscript h refers to the homogeneous case. As 
expected, it is a very small energy barrier (as compared with 
the one typically found in some optical chains, i.e., chains 
with on-site nonlinear potential, see [|20|]); for example, for a 
very heavy DB, AE^, N /E^ dd ~ 2.1%only! This explains the 
well-known fact that it is rather easy to create mobile DBs in 
an FPU chain, and also why in the first-order approximation 
in this barrier was found to be zero. In Fig. || we represent 
the dependence of the barrier on various parameters: (i) DB's 
period Tub (as expected, also see below the discussion on 
the generation of moving DBs, the PNB is larger for higher- 
frequency DBs; in the first order of the perturbative expansion, 
PNB varies as 1/T^ B ). (ii) a and (iii) /3. At the first order in 
the perturbational expansion, the PNB does not depend on a, 
but only on 1/(3, i.e., it decreases with increasing nonlinear- 
ity. This feature can be easily understood if one views the role 
of the nonlinearity as reducing particles' excursions around 
equilibrium, and therefore reducing the differences between 
the odd- and even-parity configurations, i.e., the PNB. (iv) Fi- 
nally, on m (note that, in the first order of the perturbative 
expansion, the PNB varies as to 2 ). 



B. DB generation and initialization 

For simulation purposes, the static DBs were generated nu- 
merically in the homogeneous FPU chain using the extremely 
fast algebraic method recently introduced by Tsironis J2^ ] . As 
shown in [22|, this method, although approximate, is gener- 



ally more accurate than the RWA and agrees with the exact re- 
sults of the anticontinuous limit method (which requires much 



longer computational times, see [|19|]) typically to 1% or even 
better. 

In order to move these breathers, we used a simple approx- 
imation of the systematic pinning mode excitation method of 
Chen et al. [0]. Namely, we 'kick' initially the DB by assign- 
ing to the points of the lattice initial relative velocities that are 
a fraction A (the "kicking coefficient") of the gradient of |r n |, 



i.e., 



l (t = 0) = A(|r n+ i|-|r n _i|)/2 



(14) 



Note that this method is not so different from the 'more em- 



pirical' methods used in [16| to obtain moving DBs. We no- 



tice that, when starting to move, the DB first loses, through 
phonon radiation in the lattice, a large part of the kinetic en- 
ergy we assigned to it when kicking. The rest of the received 
energy is used to overcome the Peierls-Nabarro barrier, and, 
as already mentioned, the DB moves from one lattice site to 
another by a continuous alternation between odd- and even- 
type configurations. 

This alternation between the two types of configurations for 
a moving DB can be noticed when inspecting the temporal 
evolution of the potential (or kinetic) energy of the DB. In- 
deed, the envelope of the temporal oscillations of DB's po- 
tential (kinetic) energy presents a series of periodically alter- 
nating relative maxima and minima, indicating the alternation 
between these configurations. The period between two such 
successive maxima (or minima) gives a rough estimate of the 
time needed by the DB to move from one site to another. But 
no more than a "rough estimate", because (i) the real time a 
DB takes for this movement does not bear a commensurabil- 
ity relation with Tub, and (ii) the structure of the envelope is 
more complex, due to the presence of other "secondary" fre- 
quencies of the DB (see [Q], and the discussion in Sec. V), 
and (iii) there are some "imperfections" in this periodic be- 
havior of the envelope, that are connected to the existence of a 
rather irregular time dependence of the relative phases of two 
neighboring sites (already mentioned in JI^]). Probably, this 
is ultimately related to the non-exact character of a moving 
DB as a solution of the Hamiltonian lattice. Note also that 
a moving DB constantly loses energy while moving through 
the lattice, although at a very small 'dissipation rate'. For ex- 
ample, as also shown in [24|, for a DB moving in a uniform 
lattice, this energy decrease, if fitted to an exponential, corre- 
sponds to a decay rate on the order of ~ 10~ 6 /unit time. This 
rate is higher for the faster DBs. Also, as explained in de- 
tail by the same authors, the analysis of the temporal behavior 
(and, in particular, of the extremal points) of the kinetic and 
potential energy allows one to evaluate the translational en- 
ergy of a moving DB, which was found to be at most 1% of 
the total energy of the DB. Not surprisingly, this value is of 
the same order of magnitude as the Peierls-Nabarro barrier. 

Returning to the kicking method for moving a DB, we make 
several other remarks. First that, as previously noticed (see, 
e.g., [p5|), the 'light' DB's (i.e., those that are relatively not 
too far in frequency above the phonon band limit) are def- 
initely easier to move than the 'heavy' DBs (which are, by 
comparison, much more localized and therefore much more 
sensitive to the discreteness of the lattice). In terms of the ini- 
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tial kick, this means that the minimum value of the kicking co- 
efficient, A, for which one gets an essentially regular motion of 
the DB [26 1 is larger for 'heavier' DBs. Also, one notices that, 
in general, the velocity of the moving DB obtained through 
this kicking method seems first to increase with increasing 
A, and after that it reaches a certain 'saturation value', i.e., 
it does no longer increase with A, but keeps a constant value. 
This leads to a rather narrow window of the possible values 
of the DB velocities, which is somewhere around a tenth of 
the phonon's velocity (for example, for a = j3 = 1, a lattice 
constant a = 10, and for a DB of period Tub = 2.1, the 
values of the velocities belong to a window of as [0.35, 1.25]; 
note that the sound velocity corresponds to 10 in the dimen- 
sionless units). In the simulations we used a chain with the 
first and last point held fixed (i.e., fixed boundary conditions; 
this should not raise conceptual problems, as such points cor- 
respond to m — > oo). Also, we tried to avoid the interference 
between the observed phenomena and the phonons that reflect 
on these fixed edges-and for this purpose we generally used 
sufficiently long chains (so that the reflected phonons do not 
come back to the interesting central regions during the obser- 
vation period). 



III. INTERACTION OF A DB WITH A JUNCTION 

We now address the main problem in this paper. Consider 
the junction between two semi-infinite FPU chains (let us call 
them "A" and "B", respectively, with the corresponding sub- 
scripts for their characteristic parameters). We fix the param- 
eters of the A chain a a, (3 a, and rtiA- For the B chain, we 
will fix the interaction parameters identical to those of the A 
chain, 



function: 



as = ola 
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(15) 



and vary the mass of the particles, mj. Note that in all the nu- 
merical simulations we set ctA = Pa = 1, as well as tjia = 1. 

A DB is generated in the A part of the chain and is sent 
to the junction with the B part. Depending on the difference 
between the masses of the particles in the two parts of the 
chain, the DB exhibits different behaviors at the junction. 



A. The Peierls-Nabarro barrier at the junction 

In order to understand and predict the behavior of a moving 
DB at such a junction, the first step is to study the change in 
the Peierls-Nabarro barrier at the junction. Namely, to deter- 
mine what would be the equivalent of the odd- and even- type 
configurations at the junction, and what would be the corre- 
sponding difference in the configurational energy. Note that 
in the case of an inhomogeneous chain the PNB is defined as 
the difference between the global maximum and the global 
minimum of the configurational energy. We consider the case 
when S = (wig - mA)/mA is a small quantity, that we use 
as a perturbation parameter for evaluating the changes in the 
odd and even configurations at the junction. Thus, we con- 
sider a "double" perturbation expansion of particles' envelope 
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where ea is evaluated with respect to the parameters of the A 
chain, i.e., 



m A uj db 



(17) 



Corresponding to the different configurations it has to exhibit 
in order to traverse the junction, the DB encounters three new 
energy barriers (refer to the Appendix for more details). These 
are, in the order of their appearance as the DB moves through 
the junction: 
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This energy barrier corresponds to the difference between the 
energy of the odd-type configuration I in the Appendix-for 
which the site of maximum elongation is the last site in the 
part A of the chain-and the energy of the even-type configu- 
ration in the homogeneous A chain. 
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This energy barrier corresponds to the difference between the 
energy of the odd-type configuration III in the Appendix-for 
which the site of maximum elongation is now the first site 
in the part B of the chain-and the energy of the even-type 
configuration in the homogeneous A chain). The last PNB 
barrier is associated with the difference between the energy of 
the odd-type configuration in the homogeneous B chain and 
the even-type configuration in the homogeneous A chain: 
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Here AE 



HA) 

PN 



denotes the Peierls-Nabarro barrier in the ho- 



mogeneous A chain, and j refers to the junction. 

Note that for a heavy-light junction, i.e., for S < 0, these 
barriers are smaller than the PNB barrier in the homogeneous 
A part of the chain and therefore a DB that moves smoothly 
in region A will have no 'energetic difficulties' to enter region 
B. On the contrary, for a light-heavy junction, i.e., for 5 > 0, 
these barriers are larger than the PNB in part A of the chain 
and one sees that, at the dominant orders in ea and 5, they 
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increase in succession. Therefore, there appears the possibil- 
ity that a DB that arrives at such a junction cannot overcome 
either the first, or the second , or the third barrier. The pres- 
ence of these barriers is confirmed by numerical simulations 
(through fine tuning of mj). 



B. The light-heavy junction 

We first present the generic results of our simulations for 
this case. 

(a) A DB can continue its movement into region B. Its fre- 
quency is not (detectably) modified. The DB keeps on losing 
energy in region A as well as in region B, but at a smaller rate 
in region B, see Fig. ||. This might be connected to the fact 
that (given that it keeps essentially the same frequency) the 
DB is further away from the phonon band limit in region B 
than in region A. Also, its velocity in region B is smaller than 
in region A. This is related to the fact that a part of the "ex- 
tra" energy that in region A corresponded to its movement as 
a whole (with a velocity va) is now used for the new, higher 
mean configurational energy, and also to overcome the corre- 
spondingly higher Peierls-Nabarro barrier in region B. There- 
fore, the "extra" kinetic energy, and correspondingly the ve- 
locity vb in region B are smaller than in region A. 

(b) The DB can reflect at the junction and return to the re- 
gion A. Its frequency and energy (see Fig. ||) are not sensi- 
tively modified by this reflection, and neither its velocity (that 
only changes sign). 

These observations can be explained qualitatively on the 
basis of the results presented above for the PNB that a DB 
(that keeps a constant period Tub) has to overcome in order 
to continue its movement in region B. The main conclusion 
is that for a DB arriving at the junction, there exists a critical 



value of the mass raj 



above which the DB cannot 



penetrate in region B and is reflected to region A. This critical 
value depends on the frequency of the DB, namely it increases 
with decreasing Tub (i-e-, it is larger for heavier DBs). How- 
ever, it also depends on the velocity the DB has in region A: it 
increases with increasing va (i.e., a more rapid breather needs 
a larger mass m J crit to be reflected than a slower DB of the 
same frequency). This can be readily understood: a more 
rapid DB in region A has more "extra" energy (above the 
Peierls-Nabarro barrier AEp^) than a slower one. There- 
fore, it may use this energy to overcome the Peierls-Nabarro 
barrier at the junction and to penetrate in region B, while a 
slower DB cannot overcome the junction barrier. 

We present below two comparative sets of pictures of the 
cases when a given DB (a) moves in a homogeneous chain, (b) 
passes through a light-heavy junction, and (c) is reflected at 
the light-heavy junction. Fig. fi] shows the temporal evolution 
of the configurations of the DB in these three situations, while 
Fig. ^ shows the movement of DB's center along the chain, 
and also the temporal evolution of the elongations of various 
particles affected by the DB. 



C. The heavy-light junction 

As already mentioned above, given that the PNB decreases 
at a junction with tub < m A , one can naively predict that the 
DB will always penetrate and continue to move in region B 
without any hinderance. Numerical simulations show that this 
is indeed the case-at least as long as the difference between 
tub and rriA is sufficiently small. For example, in the partic- 
ular case ofrriB = 0.99 (recall that in simulations we took 
nriA = 1) we studied the dependence of the characteristics of 
the 'transmitted' DB on those of the 'incident' one. First of 
all, one notices that the transmitted breather takes some time 
to "adjust" to the new environment (and 'heavier' DBs take 
a longer time to adjust than the 'lighter' ones). During this 
period, the DB loses energy and adjusts its final energy to the 
smaller mass of region B. 

The transmitted DB (within estimated errors) has the same 
period as the incident one (i.e., the adjustment is such that it 
preserves DB's frequency). After this transient period, the DB 
reaches a constant 'asymptotic' velocity. In general, there is 
no simple relationship between this asymptotic velocity vb 
and the characteristics of the incident DB, namely its initial 
velocity (in region A) va and its period Tub- 

However, there is a tendency towards 'uniform' veloci- 
ties after transmission through the junction for a given DB. 
Namely, for a DB with period Tub and different velocities 
in region A, va, the effect of entering region B is to reduce 
the dispersion of these velocities, i.e., the dispersion of the 
asymptotic velocities, Vb, is smaller [p7|]. This observation is 
illustrated in Fig. ^ for a given DB (with Tub = 2.1) and for 
three representative initial velocities va (chosen, respectively, 
as the lower and upper limits of the velocities that could be ob- 
tained through the "kicking" method described above in Sec. 
[IB, and one value in-between these limits). 



It was relatively more difficult to investigate the depen- 
dence of the asymptotic velocity vb on DB's period Tub, 
simply because it is rather difficult with the kicking method 
to obtain the same velocity for DBs of different frequencies. 
However, we managed to obtain four DBs of periods vary- 
ing between 2.2 and 2.5 (with a step 0.1) and almost (within 
4%) the same initial velocity. The simulations show no simple 
monotonic dependence of Vb on Tub- 

Next, we focus on the most important part of this section 
(that will clarify the meaning of sufficiently small difference 
between tub and wia)- Specifically, how does the behavior 
of a given DB depend on the value of 7713? To analyze this, 
we ran systematic simulations for a given DB (we chose one 
with Tub = 2.1 and an initial velocity va = 0.928, that 
corresponds to a kicking coefficient A = 0.7), and for vari- 
ous values of tub- The behavior of the DB at the junction is 
rather complex, and can be described essentially as follows: 
the DB, entering the region of lower mass, has "extra" energy. 
During an "adjusting period" (that might take from about ten 
to hundred DB periods), this extra energy is redistributed be- 
tween: (i) the kinetic energy of DB's translation as a whole 
(the DB is accelerated upon entering region B); (ii) perturba- 
tions (which we address later) in the A and also in the B part 
of the chain; and (iii) a slight decrease of transmitted DB's 
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period (i.e., increase of its configurational energy). The redis- 
tribution of energy between these elements is a delicate pro- 
cess, and it depends on the mass difference between regions 
A and B. When the mass difference is sufficiently small, up 
to, say, (rriA — mB)/mA = 0.4, the predominant phenomena 
are (i) and (ii)-the perturbations being small-amplitude ones, 
i.e., phonons that move rapidly far away from the junction, in 
both A and B part. 

When the mass difference is even larger, we find that in A 
part there are not simply phonons that appear, but a reflected 
DB: the initial DB, arriving at the junction, is split into a re- 
flected DB and a transmitted one. Moreover, the transmitted 
DB is usually (nonlinearly) unstable and subsequently splits 
into two (or sometimes more) other DBs. See Fig. ||for a re- 
alization of these phenomena: the trajectory of the initial DB, 
the reflected one, the transmitted DB and its subsequent split- 
ting into two other DBs. Figure ^ offers the energy variation 
associated with these phenomena. We note that the total en- 
ergies of the resulting DBs never sum up to the initial energy 
due to the phonon losses in the chain that accompany all these 
processes. To our knowledge, such DB splitting has not been 
noticed before. If we continue to decrease tub, the reflected 
DB (that is initially very weak energetically as compared to 
the transmitted one) becomes progressively more energetic, 
while the transmitted DB becomes progressively weaker and 
finally dissappears in region B leaving only rapidly-moving 
phonons in its wake. The end product is the "strong" (i.e., 
large-amplitude) reflected DB in region A. 



IV. INTERACTION OF A DB WITH A MASS RAMP' 

A. The Peierls-Nabarro barrier for a 'ramp' 

Consider that in region B the mass of the particles varies 
slightly, linearly, as one moves away from the junction point, 
i.e., the mass of the k-th particle in B part is 



rriB(k) = ?tm(1 + fcA) . 



(21) 



where A > corresponds to an up mass-ramp, while A < 
to a down mass-ramp, and for analytic calculation purposes 
we consider that |A| -c 1. A double analytical expansion 
in ea and A allows us to estimate the shape function for 
the equivalents of the odd- and even-type configurations, and 
therefore to estimate the Peierls-Nabarro barrier the DB must 
overcome in order to move up to site k in region B. The barrier 
(with details given in the Appendix) is found to be: 



AE r PN (k) = E r odd (k) - E h e ^ = AE 



HA) 
PN 



+m A uj 2 DB ( {A [(0.77906fc + 0.77906) e^ 1 

-(0.16470A; + 0.16470)] 

+A 2 [(0.38953/c 2 + 0.77906fc + 2.56207) e^ 1 

+0.91605 - 3.2585^]} , (22) 

where the superscript r refers to the ramp. 



B. The 'up-ramp' 

This corresponds to the case A > 0. The main result is 
that a DB that enters the B part of the chain is finally reflected 
(at some point within the B chain) and returns to part A. Note 
that: 

(a) The point where the DB is reflected, i.e., the critical 
mass on the ramp, m r crit , depends on the 'slope' A of the 
ramp and is generally different from the value mj,, jf tha t cor- 



responds to the case of an abrupt junction, see Sec. [II B This 
can be seen by equating the (critical values of the) most en- 
ergetic odd-type configurations in B chain in the case of an 
abrupt junction and of a ramp, and finding the relationship 
J onri m r i n me particular case shown in 
we note that the critical mass decreases with de- 



between m J crit and m r crit 



Fig. 3(a) 



creasing slope of the ramp and that it is smaller than the value 
for the junction case. 

(b) For a given ramp, the critical mass increases with in- 
creasing initial v elocit y of a given DB (with a fixed fre- 
quency), see Fig. 8(b) Sometimes the DB can get trapped, 



as seen in the inset of this figure. Note, however, that if one 
changes slightly DB's initial position in region A (without 
changing its initial velocity), then the DB is no longer trapped, 
but reflected, see the inset. Thus, trapping seems to be a rather 
delicate phenomenon, that depends on 'how' (i.e., with what 
precise configuration and relative phase difference between 
sites) the DB arrives at the trapping site. 

(c) For a given ramp, the critical mass seems to increase 
with a decrease in DB's period (i.e., it is lar ger fo r 'heavier' 
DBs for the same initial velocity), see Fig. |8(c)| . Note that 
in all these cases there is a typical temporal evolution of the 
energy of the DB. Before entering region B, one recognizes 
the usual small energy loss in an uniform FPU chain; then in 
the ramp part there is a somewhat smaller energy loss (pre- 
sumably the DB is a little bit further away from the phonon 
band) that becomes progressively smaller when the DB is de- 
celerated on the ramp. At a certain moment, the DB starts 
to 'descend' the ramp, to increase its velocity, and its energy 
loss increases progressively, again up to the usual loss in the 
homogeneous chain. When the DB 'ascends' the ramp, its 
configurational energy averaged over a period (and the corre- 
sponding Peierls-Nabarro barier) increase at the expense of its 
translational energy. Therefore, at a certain moment, the DB 
does no longer have sufficient "extra energy" to overcome the 
barrier and it is reflected (and sometimes it may get trapped). 
Rolling down the hill, it recuperates its translational energy 
and when it gets out from region B and re-enters region A it 
has almost the same velocity as its initial one in region A. The 
transmission and reflection are "almost elastic", in fact the DB 
loses a little bit less energy than it loses normally during its 
movement in a uniform chain. 



C. The 'down-ramp' 

Consider now that in the region B the mass of the particles 
decreases from one particle to another with the small quantity 
rriA A < 0. An illustration of DB's typical behavior is given 
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in Fig. ||. As the DB enters the ramp it accelerates with a 
concomitant narrowing of its shape and the emission of some 
radiation (phonons). This is clearly seen as a change in slope 
(left panel). At later times we observe another change in slope 
signifying further acceleration of the DB with significant ra- 
diation and emission of smaller breathers (right panel). One 
can also observe other secondary, small-energy DBs that may 
form at later stages. Nonetheless, these phenomena are highly 
complex and beyond our current level of understanding. 



V. CONCLUSIONS AND PERSPECTIVES 
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We have systematically explored the transport properties of 
a discrete breather in a nonlinear chain comprising two seg- 
ments with differing mass, specifically in an FPU chain. We 
considered abrupt junctions (light-heavy and heavy-light) as 
well as (up and down) mass ramps. We studied the trapping, 
reflection, transmission and splitting of the DB as a function 
of junction type, mass difference, breather frequency and ve- 
locity. The DB splitting, trapping and reflection may take 
place either at the junction or at a particular particle within 
the ramp. We also estimated the Peierls-Nabarro barrier for 
the different cases to understand the DB transport across a 
junction or within a ramp. However, the approach for calculat- 
ing the PNB is based on the fundamental assumption that the 
period (frequency) of the DB does not change "significantly" 
during its movement, which has its limitations, as shown by 
the simulations and indicated above in various cases. There- 
fore, we can rely on this method only at a qualitative level. 

In the present paper we exclusively focused on two seg- 
ments with slightly different masses. It would be interesting 
to explore a junction (or ramp) between two segments with 
the same mass but with differing strength of either the har- 
monic (a) or anharmonic (J3) interaction parameter of the FPU 
chain. This is under investigation and our preliminary results 
do not demonstrate a qualitatively different picture compared 
to the mass case. In addition, if we consider an A-B-A mass 
sandwich structure then there is a distinct possibility that the 
breather will get trapped inside the B segment. By a suitable 
choice of the mass profile one may envision a 'breather lens'. 
This is currently explored and preliminary results agree with 
these conjectures. We believe that our results are not specific 
to the FPU chain. Other nonlinear potentials should lead to 
generically similar results. Many open questions remain, e.g. 
better estimates for site-to-site traversal time of a DB, influ- 
ence of the "secondary" frequencies of the DB on its behavior 
(for example, on the envelope of temporal oscillations of en- 
ergy), a better understanding of the nonlinear instability that 
leads to DB splitting (reflected/transmitted, and afterwards to 
the secondary splitting of the transmitted DB), and conse- 
quently, to the complex behavior on a down-ramp, etc. An 
experimental realization of our findings in low-dimensional 
electron-phonon coupled materials pi2|], e.g. conjugated poly- 
mers and metal-halogen chains [H], using different iso- 
topes would be quite instructive in unraveling the interesting 
transport properties of breathers with potential applications. 



VI. APPENDIX: PNB FOR VARIOUS CONFIGURATIONS 



In this appendix we present the relevant details of esti- 
mating the Peierls-Nabarro barrier for the different cases dis- 
cussed in the text. 



A. The homogeneous chain 

1. The odd-type mode 

It is characterized by £„ = £_„ = (— 1 ) ' n, 1 77^ (rj n , the re- 
duced shape function, being positive for all n), together with 
the condition Co = 770 = 1 (that gives the normalization of the 
shape function). The equations for the reduced shape function 
read, respectively: 



2 -e(l + 77i) . 
Ae = — jt- 1 , for n = , 



(l + ^i) 3 

+ 

+ (r?„ + Vn-i) 3 } Jot n>l 



e Ae 
-(2r] n + ?7„+i + T} n _ x ) + — [(??« + ?7™+i) 3 



(23) 



Here 



A 



4a 



(24) 



Note the singular behavior in of the square of the am- 
plitude, A 2 odd (e.g., this means that the 'heavier' the DB, the 
larger its amplitude). Using the series expansions in e for £„, 
Eq. (|l2|), and the corresponding ones for the reduced shape 
function in Eq. (23), and ordering the corresponding powers 
of e, one can show that the series rp! (for a fixed j, i.e., for 
a fixed order in the perturbative expansion in e) rapidly con- 
verges to zero with increasing n; more rapidly for the case of 
small js than for larger ones [p8|. Finally, one is led to the 
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following expressions for particles' shape function: 



and the equation for the amplitude: 



£o = l, 

£1 =£_i = (-l)[0.52305 + 0.15113e + 0.08549 e 2 
+0(s 3 )] , 

£2 = £-2 = 0.02305 + 0.15691 e + 0.12643 e 2 

+0(e 3 ), 

£ 3 =£_ 3 = (-1) [0.00580 £ + 0.04239 e 2 
+0(e 3 )}, 

f„ = 0(HT 6 ) , |n| > 4 , (25) 

together with the dependence of the amplitude A dd on DB's 
frequency, mass m of the particles (through e), a, and f3: 

A = H ' odd = 0.56609E" 1 - 0.59960 + 0.02366 e+0(e 2 ). 

(26) 

All these lead finally to the following expression for the con- 
figurational energy of the odd-parity mode: 



E odd = mUJ 



DB 



(3 



[0.38953E- 1 -0.16470 



0.12386 e+0{e 2 )} 



(27) 



2. The even-type mode 

It is characterized by the presence of two 'main peaks', 
£0 = — £1 = 1, and also by a staggered shape: £„ = 
— £-n+i = ( — l)'"'>7n> w i m m e reduced positive shape func- 
tion rj n . In this case, the equations for the reduced shape func- 
tion read, respectively: 

Ag= 4- 1 (3 + m ) 



8+(l + r?i) 3 
+(r?« + ?/«-i ) 3 ] ,for n > 2. 



£ As 
??„ = — (2r7„ + ?7„ + i + r/„_i) + — [{Vn + Vn+i) 3 



Here 



^ 3/?A et;en 

4a 



(28) 



(29) 



Finally, one is led to the following expressions for particles' 
shape function: 

Co = -6 = 1 , 

6 = = 0.16579 + 0.31767e + 0.13806 e 2 + 0{e 3 ) , 
£3 = -£_ 2 = (-1) [0.00048 + 0.04438 £ + 0. 10766 e 2 

+0(e 3 )] , 
£ 4 = -C-3 = 0.00012 e + 0.01115 e 2 

+0(e 3 ) , 

£ n = 0(HT 6 ) ,n>5, n<-4, (30) 



A = 



3/M 2 „ e?I _ 0.41735 
4a £ 



0.38670 + 0.02077e + C(e 2 ) . 

(31) 



The corresponding configurational energy: 



El 



m ^DB ( 7; ) [Q.38117^ 1 - 0.15705 



- 0.10559 £ + C(£ 2 )]. (32) 

Note that E% ven < E^ dd , i.e., as already noticed [|l^|, the 
even-type mode is more stable than the odd-type one. 



B. The junction 

We refer to Fig. [l(] to follow the different configurations 
of the DB moving from left to right through the junction. As 
indicated in the text, Eq. (|l6|), we used a double perturba- 
tion expansion of the envelope function-in both ea (evaluated 
with respect to the parameters of the A chain, see the text) and 
8 to compute the different configurations. Note that the con- 
vergence in 8 is not as good as that for ea\ namely, 8 should 
be 10 2 times less than ea in order to get the same degree of 
correction for the same order of expansion in S as in e A . How- 
ever, the details of the calculations are lenghty and, because 
they present no conceptual difficulty, not given here. Instead, 
we give the expressions for the configurational energies of the 
DB in its successive appearances-as these allow us to com- 
pute the various Peierls-Nabarro barriers it encounters. 



I. Configuration I 

It is of the odd-type-it corresponds to the first panel in 
Fig. [Kj: the site of maximum elongation is in part A of the 
chain. Its energy is found to be: 



E 



odd 



E 



h(A) 
odd 



m A u DB 



OA 

Pa 



x [8 (0.13793 e~ a x - 0.01573 + 0.0056 e A ) 

+8 2 (0.49933 e~ a x + 0.21657 - 0.69346 e A )] . (33) 

Correspondingly, the first barrier that the DB has to overcome 
is the one between an even-type configuration in the homoge- 
neous A chain and this configuration, namely: 



E 



m 

odd 



p h(A) _ a p h{A) , 2 A 

^everi ~ ^ E PN + m AU DB \ — 



x [8 (0.13793 £4* - 0.01573 + 0.0056 £,4) 

+8 2 (0.49933 e a 1 + 0.21657 - 0.69346 e A )] . (34) 

2. Configuration II 

It is of an even-type and corresponds to the second panel of 
Fig. [HI there are two sites with large elongations, one of them 
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in part A of the chain, the other one in part B of the chain. The 
corresponding energy, up to 0{e 2 A , 8 2 ), is: 



and, on the other hand, the PNB in the homogeneous B part 
is: 



a a 
Pa 



x [5 (fJ.38117^ 1 - 0.07852) 

+<5 2 (-0.83444 e^ 1 - 0.32507+ 1. 11672 e A )] ■ (35) 

Therefore, the energy difference between configurations II 
and I is: 



E h(B) _ Tjjh(B) _ A E HA) 



2 I a A 
n A^ DB [ -j^ 



x [(5(0.01672 e^ 1 - 0.00765) + S 2 (0.00836 e^ 1 )} . (42) 



C. The ramp 



FA") - K j{1) - -AF h{A) 1 - 2 1 aA 



odd - ~ LJl ^PN "T" m -A UJ DB y-p^ 

x [5 (-0.24324 s A x - 0.06280 - 0.00526 e A ) 

+S 2 (-1.33278 e^ 1 - 0.54163 + 1.81108 e A )] .(36) 



Consider that the ramp has a "slope" A, i.e., the mass of the 
k-th site in the ramp is 



ms(fc) = m A (l + kA) 



(43) 



3. Configuration III 

Again of the odd-type, it corresponds to the third panel in 
Fig. [lCj the site of maximum elongation is now in region B of 
the chain. Its energy: 



^odd ~ 



x [S (0.64113 e A x - 0.14880 - 0.00526 e A ) 

+S 2 (0.75093 e A x + 0.21657 - 0.68820 e A )} . (37) 

Therefore, the energy barrier between configuration II and 
configuration III is: 



E 



idd 



EiVi) - AE h{A) 

^ even ~ ^-^PN 



2 / a A 
A^DB [ 



: [S (0.25996 i 
-5 2 (1.58537 e; 



A 



0.07045- 0.00526 e A ) 

0.54163- 1.80493 e A )} .(38) 

After this, the DB is essentially in the homogeneous B part 
[of mass tub = + 6)]; the energies of the even- and 

odd-type configurations in B are: 



F h(B) F h{ A ) 



m A LO DB 



OA 

Pa 



[5 (0.76234 e A x - 0.15705) 
+S 2 (0.38117 e^ 1 )] , 



h odd - h odd + m A UJ DB \-p^ 

x [5 (0.77906 e A x - 0.16470) 
+5 2 (0.38953 e^ 1 )] . 



(39) 



(40) 



Therefore, on one side, the energy difference between the con- 
figuration III and the even configuration of chain B is: 



F KB) _ p. 
^even odd 



3(H) 



-AE 



HA) 
PN 



m A uj DB 



OA 

Pa 



[S (0.12122 e^ 1 - 0.00808 + 0.00526 e A ) 
-S 2 (-0.36976 e~ A l - 0.21657 + 0.68820 e A )\ ,(41) 



Note that A > corresponds to an up-ramp, while A < to 
a down-ramp. Consider first a configuration of the even-type 
where the two sites with maximum elongation are k — 1 and 
k. The corresponding configurational energy is found to be: 



e: 



E, 



KA) 



2 / a A 
A^DB ( -gT" 



x {A [(0.76234fc - 0.38117) e^ 1 
+ (-0.15705fc + 0.07852)] 
+A 2 [(0.38117fc 2 - 0.38117fc - 1.35992) e 
-0.85452 + 3.28732 e A ]} . 



-l 

A. 



(44) 



Consider then the next configurational step in the displace- 
ment of the DB from left to right on the ramp, i.e., an odd- 
parity type of configuration that is centered on the fc-th site, 
i.e., the /c-th site is the one that has the maximum elongation. 
The energy of this configuration is: 



Kddik) 



E 



KA) 

>dd 



2 / a A 
A^DB I 



x {A [(0.77906fc + 0.77906) e~ A 

-(0.16470fc + 0.16470)] 

+A 2 [(0.38953fc 2 + 0.77906/c + 2.56207) e~ A 

+0.91605 - 3.25850 e A ]} . (45) 



Correspondingly, the energy barrier to overcome while 
moving from site k — 1 to site k is: 



E odd(k) - Elven(k) 



AE 



HA) 
PN 



2 / a A 

A DB ' ~Ba 



x {A [(0.01672fc + 1.16023) e~ A 

-(0.00765/c + 0.24322)] 

+A 2 [(0.00836A; 2 + 1.16023fc + 3.92199)e^ 1 

+1.77057- 6.54584 e A ]} . (46) 

At the dominant order in both A and e A , one finds an increase 
in the PNB for an up mass-ramp (A > 0), and a decrease of 
the barrier for a down mass-ramp (A < 0). 
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behavior of the coefficients of the perturbative series intuitively 
more accessible. 



n-4 n-3 n-2 n-1 11 n+1 n+2 n+3 n+4 n+5 
site number 



(a) homogeneous chain 




(b) light-heavy junction, ms = 1.002: transmission of the DB 
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(c) light-heavy junction, niB = 1.04: reflection of the DB 



FIG. 1: Temporal evolution of a DB's configuration in three differ- 
ent situations. Notice the alternation between odd and even- type 
configurations. Tub =2.1. 



FIG. 2: Variation of the Peierls-Nabarro barrier as a function of 
breather time period (Tdb), FPU chain parameters (a, j3) and mass 
(m). (a) a = 1, = 1, m = 1. (b) T DS = 1.2, /3 = 1, m = 1. (c) 
a = 1, T DB = 1.2, m = 1. (d) a = 1, p = 1, T DS = 1.2. 
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(a) homogeneous chain 





(c) light-heavy junction, rriB = 1.04: reflection of the DB 



FIG. 4: Propagation, transmission and reflection of the DB for the 
three situations in Fig. |l|. The associated motion of particles in the 
region of the DB is also shown. Tdb = 2.1. 
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FIG. 5: A heavy-light (tub = 0.99) junction. For various initial 
velocities da of a DB with a given period (Tdb = 2.1) one gets 
approximately the same asymptotic velocity vb in region B (see the 
lower inset). The upper inset depicts the DB's position at very early 
time before reaching the junction, while the lower inset depicts the 
DB's asymptotic trajectory. 
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FIG. 6: A heavy-light (ms = 0.50) junction. An initial DB 
{Tdb = 2.1) splits into a reflected and a transmitted DB. Later on, 
the transmitted DB furhter splits into two other DBs (the circles on 
the figure indicate the regions of the splittings). Note that when the 
mass difference is even larger, the transmitted DB (that has progres- 
sively less energy) might split into three or even four smaller DBs. 
Decreasing ms furhter leads practically to the disappearance of the 
transmitted DB, and to a substantial phonon creation. 
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FIG. 7: The energy associated with the phenomena described in 
Fig. ^j. The dashed-dotted lines on the figure delimit the intervals of 
the occurrence of the splitting phenomena, when there is no net sep- 
aration between the resulting DBs, i.e., no clear separation of their 
energies. 
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FIG. 8: The behavior of a DB on an up-ramp depending on (a) 
the slope of the ramp (Tdb = 2.1) and (b) its initial velocity va- 
Note that the DB can also get trapped on the ramp; but, as shown 
in the inset, a slight perturbation-for example, a slight modification 
of the initial conditions-can lead to the disappearance of trapping. 
(Tdb = 2.1). (c) The behavior of a DB as a function of the pe- 
riod Tdb of the DB. (Note the limited possibilities to obtain DBs of 
various periods with rigorously the same initial velocity va) 
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FIG. 9: Acceleration of a DB on a down-ramp. Site n=0 corresponds 
to the last site in the region before the onset of the ramp. The slope of 
the ramp is S — 0.0018. On the y axis we approached the sites by 9.9 
lattice constants (a = 10 units) in order to increase the resolution. 
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FIG. 10: Schematic representation of the succession of the odd- and 
even- type of configurations for a DB traversing a junction from left 
to right. The white and black circles correspond, respectively, to 
particles in the A and B parts of the chain. 



